Draft version November 17, 2008 

Preprint typeset using 1^1^^ style cmulatcapj v. 10/09/06 



DISCOVERY OF FULMINIC ACID, HCNO, IN DARK CLOUDS^ 
NuRiA Marcelino^, Jose Cernicharo^, Belen Tercero^, and Evelyne Roueff^ 

Draft version November 17, 2008 

ABSTRACT 

We report on the first detection in space of fulminic acid, HCNO. This isomer of HNCO has been 
observed in three starless cores, Bl, L1544 and L183, and in the low mass star forming region L1527 
with a measured abundance ratio of HNCO/HCNO between 40-70. However, HCNO was not detected 
towards the direction of the cyanopolyyne peak of TMC-1 or towards the Orion Hot Core region. The 
derived HNCO/HCNO abundance ratio in these cases is greater than 350 and 1000 in TMC-1 and 
Orion, respectively. We find that CH2 -I- NO HCNO -f H is a key reaction for the formation of 
ful minic acid. A value of 5.5xl0~^^ cm^s~^ of the corresponding reaction rate coefficient, as given 
bv iMiller et al. I ()2003l ). allows to reproduce the observed abundances of fulminic acid in both the 
observed dark clouds and low mass star forming core, where the determined abundance of HNCO in 
these regions with respect to molecular hydrogen is l-5xl0~^°. 

Subject headings: astrochemistry — line: identification — ISM: abundances — ISM: clouds — ISM: 
molecules 



1. INTRODUCTION 

Isocyani c acid (HNCO) was fi rst detected towards 
Sgr B2 (jSnvder fc Buhl I Il972f l and its rotational 
transitions have been found to be prominent to- 
wards the molec ular clou ds of the Galactic Center 
(iKuan fc Snvder I Il9 96: C ummins et al. I Il986t iTurner I 
II99II : iMartm et al. I I2008f ). Since then, it has been 
observed towar ds dense cores assoc i ated with massive 
star fo rmation (iJackson et al. 1 119841 : IChurchwell et al. 1 
119861 : IZinchenkg et al. I |2000D. in the direction of 
TMC-1 (iBrownl [l98l[ iKaifu et all [200l. in diffuse 
clouds ([Turner et al. I 119991) and in exter nal galaxies 
(jMeier fc Turner II2005I : iMartfn et al. II2006I ). Hence, iso- 
cyanic acid is found in a large variety of physical envi- 
ronments. In the 3 mm line survey we have performed 
in four dark cloud cores using the IRAM 30-m telescope 
(jMarcelino et al. Il2007al lbl) we detected HNCO in all the 
observed sources, showing that this species is a common 
constituent of dark clouds. While its abun dance in dark 
and diffuse clouds is ~10~^° (jTurner Ill991h . in hot cores, 
the Galactic Cent er clouds and in externa,l galaxies it is 
~ 1-5 X 10"^ (IZinchenko et al. 1 12OOOI : iMartfn et all 
[20M [2OO8I : [Meier fc Turner II2005D . 

HNCO has several isomers with a singlet electronic 
state: HOCN, HCNO, HONG, and HNOC. The chemi- 
cal pathways leading to those isomers could be very dif- 
ferent, since there is not an obvious common precursor. 
Thus, their detection in space could provide some infor- 
mation on the chemistry of these species. Unfortunately, 
only HNCO and HCNO have been fully characterized in 
the spectroscopic laboratories so far (see below). In this 
Letter we report on the first detection of fulminic acid, 
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HCNO, towards Bl, L1544, L183 and L1527. 

2. OBSERVATIONS AND RESULTS 

The observations were performed using the IRAM 30- 
m telescope (Granada, Spain) between 2002 December 
and 2007 February, and in 2008 July. We used two 
3 mm SIS receivers working simultaneously at the same 
frequency, but with orthogonal polarizations. Both re- 
ceivers were tuned in single sideband mode with image 
rejections '^24 dB. System temperatures were ~150 K, 
for the lowest frequencies and between 200-250 K for 
the highest ones. The observations were done in fre- 
quency switching mode. The spectrometer was an au- 
tocorrelator with 40 kHz of spectral resolution (^0.13 
kms~^). Intensity calibration was performed using two 
absorbers at different temperatures. The atmospheric 
opacity was obtained from the measureme nt of the sky 
emissivity and the use of the ATM code ()Cernicharo I 
Il985f ). Weather conditions were typically average sum- 
mer conditions (water vapor column ~5 mm and zenith 
opacities <0.1 at 3 mm), except for the period in 2008 
July when we had very good weather conditions (<2 mm 
of precipitable water and opacities '^0.03). 

Pointing and focus were checked, on strong and nearby 
sources, every 1.5 and 3 hours respectively. At the ob- 
served frequencies the bcamwidth of the antenna is in 
the range 26" — 22" and the main beam efficiency is 
0.77 — 0.74. The sources being similar in size, or slightly 
larger, than the main beam of the telescope, the con- 
tribution of the error beam to the observed intensities 
is negligible and all the spectra have been calibrated in 
main beam temperature scale. The observed lines are 
shown in Figure [T] and the derived line parameters, ob- 
tained from gaussian fits using the GILD AS package^, 
are resumed in Table [TJ 

During the 3 mm line survey (85.9-93.1 GHz) of the 
dark clouds Bl, L1544, L183 and TMC-1, wc have de- 
tected several unidentified lines. Among them one line 
at 91.751 GHz was observed toward all sources except 

* http://www.iram.fr/IRAMFR/GILDAS 
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TABLE 1 

HNCO AND HCNO OBSERVED LINE PARAMETERS 
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Note. — Number in parentheses are Icr uncertainties in units 
of the last digits. Adopted rest frequencies (MHz) for the observed 
Unes are 91751. 320±0.004 & 114688.383±0.005 for HCNO J = 4-3 
and 5 - 4 respectively, and 87925.240±0.030 & 109905.757±0.030 
for HNCO 4o4 — 3o3 and 5o5 — 4o4 respectively. 

TMC-1. The observed frequency was not found in th e 
availabl e molecular catalogs: JPL (iPickett et al. II1998D. 
CDM S (iMhiler et al. I [200IL l2005h . and NIST (|Lo^ 
II992I I2004D . Using a personal spectral line catalog devel- 
oped by one of us (J. Cernicharo), we tentatively identi- 
fied this line as the transition J = 4 — 3 of fulminic acid 
(HCNO). The frequencies of this species have been mea- 
sured in the laboratory u p to Juv = 12, corresponding 
to i^max = 275.2 GHz, bv IWinnewisser fc Winnewisser I 
(|l971). The agreement between the frequency measured 
in the laboratory for the J = 4— 3 transition and that ob- 
served in dark clouds was better than 30 kHz, i.e., similar 
to the uncertainty of the laboratory measurements and 
of the velocity determination in these objects. Within 
the survey frequency range we have also observed the 
4o4 — 3o3 transition of HNCO in all sources, including 
TMC-1. Frequencies for HNCO have been computed 
fro m the effective r otatio nal constants obtained recently 
bv lLapinov et al. I (|2007t ). 

In order to confirm the new detection we had additional 
observations with the IRAM 30-m radio telescope in 2008 
July. We improved the S /N ratio on the J = 4 — 3 tran- 
sition, and observed the J = 5 — 4 towards the sources 
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Fig. 1. — Line profiles for the observed transitions of HNCO (left 
panel), and HCNO (right panel). Grey lines represent the HNCO 
(4o4 — 3o3) and HCNO (4 — 3) lines, while the black ones are for 
the HNCO (5o5 - 4o4) and HCNO (5 - 4) transitions. In Bl, the 
hyperfine components of HNCO (4o4 — 3o3) at 4 and 9 kms~^ are 
shown. Note the non detection of HCNO (4 - 3) towards TMC-1, 
and the presence of a unidentified line in this source, L1544 and 
Bl at 91750.68±0.03 MHz. 



where the previous one was detected. We also included 
in this run L1527 to our targets. This source hosts a 
Class O/I protostar and lately has been of great inter- 
est because of its strong lines of carbon-chain species 
(jSakai et al. I l2007l I2008D . In this source we have ob- 
served two positions, the central one towards the direc- 
tion of the protostar, and another one located at an off- 
set of (20", —20" ), where the emissio n of carbon chains is 
still intense (see. lSakai et al. Il2008f ). The observation of 
these two positions in L1527 could permit to distinguish 
between the emission in the region around the protostar 
and that of the cold envelope. We have also observed 
in this run another transition arising from HNCO, in or- 
der to derive volume densities, column densities, and the 
abundance ratio between both isomers. 
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TABLE 2 

Derived column densities for HCNO and HNCO. 



Source 




iV (^rllNOUj 




-ft 




IQio cm-2 


1012 cm-2 


10^ cm~'^ 




Bl 


21(1)/17(4) 


9{l)/8(2) 


6(1) 


42/47 


L1527 


6(1)/4{1) 


2.2(6)/1.8{4) 


5(1) 


37/40 


L1527-B 


5(1)/4{1) 


0.9(5)/0.8{2) 


5(1) 


18/20 


L1544 


8(3)/6{2) 


5{1)/4(1) 


0.7(2) 


62/66 


L183 


6(2)/5{l) 


3.4{3)/2.5(6) 


0.5(1) 


57/50 


TMC-1 


<1.3/<1.4 


5.7(4)/4(l) 


0.3(1) 


>390/320 


Orion-KL 


<890 


9000(2000) 




>1000 



Note. — _R is the HNCO/HCNO abundance ratio. Errors to 
the column density are in parentheses in units of the last digit. 
L1527-B corresponds to the position (20",— 20") of this source. 
The first entry in Af(HNCO), Af(HCNO), and R corresponds to 
the values obtained from the rotational diagrams, while the second 
one corresponds to those obtained from the LVG calculations. 

Figure [T] shows the observed profiles for the HNCO and 
HCNO transitions. The J = 4 - 3 transition of HCNO 
is clearly detected in Bl, L1544, L183 and L1527 weU 
above 6cr level. However, it was not detected towards 
the direction of the cyanopolyyne peak of TMC-1, with 
an rms of 6 mK. The J = 5 — 4 transition was detected 
only towards Bl and L1527 with similar intensities than 
those of the 4 — 3 line. Upper limits (3cr) to the J = 5 — 4 
line intensity are <0.010 K, <0.012 K, and <0.018 K for 
L1527 (20", -20"), L1544, and L183 respectively. This 
line should be 2-2.5 times weaker than the J = 4 — 3 
transition, which is consistent with the intensity ratios 
for HNCO in these sources. 

3. DISCUSSION 

The dipole m oment of fulminic acid is 3.099 D 
(|Takashi et al.. I [1989 ,) . while that of HNCO is 2.1 
(Aia=1.60 D, /ib=1.35 D; iHocking et a"nil97l . Con- 
sequently, both molecules could have different excita- 
tion conditions and we have performed rotational dia- 
grams in order to obtain rotational temperatures (Tj-ot) 
and column densities [N). Two transitions of HNCO 
were observed in all the sources. The fit of both lines 
gives rotational temperatures of 12±3 K, 12±4 K, 
15 ± 6 K, 6 ± 1 K, 7 ± 1 K and 6 ± 1 K for Bl, L1527 
(0,0), L1527 (20", -20"), L1544, L183 and TMC-1, re- 
spectively. Errors to the derived temperatures and col- 
umn densities include those obtained from the gaussian 
fit, calibration errors and the fact that we are using only 
two transitions to perform the rotational diagram. Al- 
though HNCO presents hyperfine structure, we have not 
included it in our calculations since it is only important 
for low rotational quantum numbers (J < 3). Never- 
theless, two weak hyperfine lines are visible towards Bl 
in the 4o4-3o3 transition at 4 and 9 kms^^ (see Fig. [1]). 
They are ten times weaker than the strongest component 
indicating that the HNCO lines are optically thin. Bl 
and the central position of LI 527 are the only sources 
where two transitions of fulminic acid have been de- 
tected, and where we have estimated rotational tempera- 
tures of 10±4 K and 17±8 K respectively. For the other 
three sources Tj-ot has been assumed to be the same for 
both isomers. Table [2] shows the obtained column den- 
sities for HCNO and HNCO, and the abundance ratio 
between both isomers. 



Collisional rates for HNCO are available from Green's 
calculations ^. We have used a Large Velocity Gra- 
dient code to derive volume and column densities for 
HNCO and HCNO in order to better constraint the 
HNCO/HCNO abundance ratio (hereafter refered to as 
R). The collisional rates J ^ J' for HCNO are un- 
known but we have assumed that they are identical to 
the Joj — > Jqj/ collisional rates of HNCO. The adopted 
kinetic temperatures are 10 K for all dark clouds except 
Bl and L1527 for which we have assumed Tkin = 15 K 
following the rotational temperatures derived above from 
HNCO. The obtained values for ^(Ha), iV(HNCO) and 
iV(HCNO) are also given in Table [1 Errors to the LVG 
calculations are estimated to be about 25 %. The agree- 
ment between the column densities derived from the two 
methods is very good. 

In Orion-KL we have detected more than 60 spectral 
lines of HNCO at 3, 2 and 1.3 mm (see the line survey of 
Orion by Tercero et al., 2008, in preparation). Column 
densities have been calculated using LTE approximation 
included in radiative transfer codes developed by J. Cer- 
nicharo (Cernicharo 2008, in preparation). We have kept 
constant the physical properties of each component of 
Orion-KL (e xtended ridge, coni pact ridge, plateau and 
hot core; see lBlake et al. II1987I for the standard values 
we have assumed) including their sizes and offsets with 
respect IRc2 (survey pointing position). Corrections for 
beam dilution are applied for each line. We have es- 
timated the uncertainty to be about 25 %, considering 
error sources as the line opacity effects, the simplifica- 
tion of the model, the division of the whole line emission 
in different components and pointing and calibration er- 
rors. The results are (1.5±0.4) x lO^^^, (1.5±0.4) x 10l^ 
(1.0±0.3)xl0^^ and (7±2)xl0i'* cm"^ ^j. ^j^g extended 
ridge, the plateau, the compact ridge and the hot core 
respectively. For the hot core a second component with 
rkin=300±75 K and 7V(HNCO)=(6±2) x lO^^ cm'^ is 
needed to reproduce the observations. In none of these 
components HCNO has been detected. The upper limit 
to R is 1100 in the hot core, the plateau and the compact 
ridge, and 300 in the extended ridge. Since the Orion 
region is dominated by grain-surface chemistry -which 
explains the large column density of HNCO- rather than 
ion-molecule and neutral- neutral gas phase reactions, the 
lack of HCNO suggests that grain surface chemistry is 
not the main path to the production of HCNO in molec- 
ular clouds (see below). 

From recent state of the art ah initio calculations, the 
relative energies of HOCN, HCNO and HONG, in their 
ground electronic singlet states, compared to HNCO, are 
24.7, 70.7 an d 84.1 kcal mol~^ at an u nprecedented level 
of accuracy (jSchuurman et al~l 120041. Another isomer, 
HNOG, is predicted b y ab initio calculatio ns at ~135 kcal 
mol~^ above HNCO (|Mebel et al. 1119960 . The relative 
abundance of the different isomers in molecular clouds 
depend on the chemical paths leading to their forma- 
tion. Unlike the case of HCN/HNC, which arc mainly 
formed from the dissociative recombination of HCNH+ , 
there is not an obvious common parent molecule leading 
to the simultaneous formation of the various HNCO iso- 
mers. The most evident formation path for HNCO and 
HCNO could be the reaction of HCNH+ with atomic 

^ unpublished results, see .http: / / molscat.giss.nasa.gov| 
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TABLE 3 

Model abundances obtained for the HNCO isomers in cold 

CLOUDS 



n(H2) 


NO 


HNCO 


HCNO 


HOCN 


NCO 


R 




10-^ 


10-10 


10-" 


10-13 


10-8 




0.1 


0.005 


0.004 


0.04 


0.47 


0.02 


0.8 


0.3 


0.01 


0.01 


0.08 


1.30 


0.07 


1.5 


1.0 


8.6 


18 


6.30 


38.0 


6.6 


29 


3.0 


7.4 


18 


3.10 


12.0 


4.5 


58 


10 


5.1 


9.3 


1.20 


2.0 


2.0 


72 


30 


3.2 


3.4 


0.58 


0.31 


0.8 


58 


50 


2.5 


2.0 


0.41 


0.13 


0.5 


49 


100 


1.7 


0.9 


0.25 


0.04 


0.3 


36 



Note. — R is the HNCO/HCNO abundance ratio 

oxygen and OH. Unfortunately, it has been shown in 
the l aboratory th at HCNH+ does not react with oxy- 
gen (jScott et al. Ifr999) an d noth ing is known about its 
reaction with OH. Iglesias] (|1977f ) suggested the reaction 
NCO+ + H2 ^ HNCO+ + H as a starting point for the 
chemistry of HNCO. The other isomers can not be pro- 
duced via this reaction as the corresponding channels are 
endothermic. We have introduced a tentative gas phase 
chemical network with the three most stable isomers - 
HNCO, HOCN and HCNO- and their as sociated ions 
which have been studied theoretically by iLuna et al. I 
()1996[ ). We have obtained that the proton transfer re- 
actions involving H;^ and NCO, the most stable isomer, 
could lead to both HCNO+ and HOCN+ via exother- 
mic channels. The HNCO+, HCNO+ and HOCN+ ions 
can further react with molecular hydrogen, giving respec- 
tively H2NCO+, H2CNO+ and HOCN H+, who s e rela - 
tive energies have been computed bv Hop et al. I (119891 ). 
These ion-molecule reactions are followed by dissocia- 
tive recombination where we assume that no chemical 
structural change occurs so that H2NCO+/H2CNO+ + 
e" HNCO/HCNO + H, and HOCNH+ + e- 
HNCO/HOCN + H. 

Neutral-neutral reactions may also be at work. So, 
HNCO can be formed by CN O2 ^ NCO + O, fol- 
lowed b y NCO + B2 HNCO -I- H, as suggested by 
[Turner I (|2000l) . However the presence of an activation 
barrier in the NCO -I- H2 reaction can not be dismissed 
and a tentative barrier of 500 K is introduced in the 
present work. In addition, a favorable neutral-neutral 
pathway leading to the formation of HCNO is provided 
by CH2 + NO ^ HCNO -I- H, a reaction we ll studied in 
comb ustion and atmospheric chemistry (.Glarborg et al. 1 
|1998|). for which no activation barrier has been found 
bv lRoggenbuck fc Tem^l (|1998D from ab initio studies. 
Other neutral-neutral reactions of HCNO in the pres- 
ence of atomic oxygen and NO may occur (|Miller et al. I 
I2OOS) . Including these various possibilities in a chem- 
ical network, we have computed equilibrium solutions 
by solving directly the steady state chemical equations. 
Table [3] displays the fractional abundances relative to 



molecular hydrogen for different volume densities in cold 
clouds, at a temperature of 10 K and a cosmic ionization 
rate of 5xl0~^^ s. These results do not differ consider- 
ably when assuming a temperature of 15 K and 20 K. The 
predicted abundance for NO is in g ood agreement with 
the v alue observed in dark clouds (jCerin et al. 1 119921 
Il993f ). NCO has also a large abundance, but due to its 
moderate dipole moment (0.67 D) and a ^H structure di- 
luting the emission spectrum resulting from an enhanced 
partition function, the expected intensity from its rota- 
tional l ines is rather weak. From our line survey of dark 
clouds (|Marcelino et al. Il2007allbl ) we derive limits to its 
abundance of l-5xl0~^. The result from our models 
show a good agreement with the values observed for R 
in dark clouds and low-mass star forming regions. It 
increases with volume density up to n(H2)=10^ cm-^, 
but decreases again when the density is larger than this 
value. It seems that for low and high volume density 
clouds fulminic acid could be rather abundant, while in 
clouds of moderate density, like TMC-1, it will be much 
less abundant than isocyanic acid, HNCO. Furthermore 
TMC-1 shows a very peculiar chemistry, being a car- 
bon rich source while it is deficient in oxygen bearing 
species. That could be the reason why HCNO is de- 
tected in L183, similar to TMC-1 in density and temper- 
ature but oxygen rich, and not towards TMC-1. The 
formatio n of HNCO on dust grains may also be very 
efBcient ()Hasewaga fc Herbst I I1993D and these surface 
chemistry processes could also lead to the formation of 
HCNO (Herbst, private communication). However, the 
gas phase models predict abundances for HNCO very 
close to those observed in dark clouds. Moreover, these 
models also predict the observed HNCO/HCNO abun- 
dance ratio. The detection of the different isomers of 
HNCO provides a fundamental information concerning 
the physical and chemical conditions, in particular the 
respective contribution of gas phase and grain processes 
in the chemistry since there is not an obvious common 
precursor for them. Thus, the observation of the other 
isomer of HNCO, HOCN, for which present gas phase 
models predict a very low abundance in dark clouds could 
provide a good discrimination between gas and grain sur- 
face processes in dark clouds. The search for HNCO iso- 
mers towards different environments would be the next 
step to further investigate their chemical production. 
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